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We study the collective excitations of a neutral atomic Bose-Einstein condensate with gravity-like 
1/r interatomic attraction induced by electromagnetic wave. Using the time-dependent variational 
approach, we derive an analytical spectrum for monopole and quadrupole mode frequencies of a 
gravity-like self-bound Bose condensed state at zero temperature. We also analyze the excitation 
frequencies of the Thomas-Fermi-gravity (TF-G) and gravity (G) regimes. Our result agrees excel- 
lently with that of Giovanazzi et al. [ Europhysics Letters , 56, 1 (2001)], which is obtained within 

' the sum-rule approach. We also consider the vortex state. We estimate the superfluid coherence 

t^r^ \ length and the critical angular frequencies to create a vortex around the 2-axis. We find that the 

. TF-G regime can exhibit the superfluid properties more prominently than the G-regime. We find 

that the monopole mode frequency of the condensate decreases due to the presence of a vortex. 

^ ; 

Ph, PACS numbers: 03.75.Fi, 04.40.-b, 34.20. Cf 

^ ■ 

^ ■ I. INTRODUCTION 

> ■ 

' The discovery of Bose-Einstein condensation (BEG) in a dilute alkali atomic gas opens up new perspective in the 
l> ! field of many body physics Most of the properties of these dilute gas can be explained by considering only 

£12 ' two-body short range interaction which is characterized by the s-wave scattering length a Recently, a new 

kind of atomic BEG has been proposed by D. O'Dell et al. [Q. They have shown that the particular configuration 
of intense off-resonant laser beams gives rise to an effective gravity- like 1/r interatomic attraction between neutral 
' atoms located well within the laser wavelength. This long range interaction potential is of the form, U{r) = —u/r, 
"Y"^ ' where u = (ll7r/15)(/aQ/cegA|). / and Al are the total laser intensity and wave length respectively, cio is the atomic 
' polarizability at the frequency 2hc/Xl. In this system, the gravity-like 1/r attraction balances the pressure due to the 
, zero point kinetic energy and the short range interaction potential. For strong induced gravity-like potential, the BEG 
I ' becomes stable even without the external trap potential |5[| . There is a competition between the gravity- like potential 
' ^ \ either with the kinetic energy (G) or the two-body short range interatomic interaction potential characterized by the 
• s-wave scattering length a (TF-G), which gives two new regimes for new atomic BEG. These two regimes are obtained 
^ ' based on the Gaussian variational ansatz for the ground state wave function Q. 

• • In the TF-G regime, collective excitation frequencies has been calculated numerically by solving the equations 

of coUisionless hydrodynamics Moreover, in this regime, an analytic expression of the ground state density is 

^ \ obtained |]. Within the sum-rule approach, collective excitation frequencies of a gravity-like self-bound Bose gas has 
' been discussed in Ref. |^]. There has been no systematic and detail study on the collective excitation frequencies 

5^ , and vortices of a gravitationally self-bound Bose gas by using the time dependent variational approach. The purpose 
of this paper is to give an analytic expressions for collective excitation frequencies, superfluid coherence length and 
critical angular frequencies required to create a vortex of a rotating Bose condensed state and to compare qualitatively 
the results of the TF-G regime with the results obtained in the TF regime of an ordinary atomic BEG. 

Here, by using the time-dependent variational method, we obtain the excitation spectrum of a gravity-like self- 
bound Bose gas. We also calculate the lower bound of the ground state energy, monopole and quadrupole mode 
frequencies of the TF-G and the G regimes. Our result agrees excellently with that of ref. which is obtained within 
the sum-rule approach. Next, we consider a rotating Bose condensate state with a single vortex along the z-axis. We 
estimate the superfluid coherence length and the critical angular frequencies required to create a vortex along the 
2-axis. We find that the TF-G regime of a gravitationally self-bound Bose condensed state should exhibit superfluid 
properties prominently than the G-regime. We find that the monopole mode frequency of the condensate decreases 
due to presence of the vortex. 
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II. TIME-DEPENDENT VARIATIONAL ANALYSIS 



The equation of motion of the condensate wave function is described by the generahzed Gross-Pitaevskii equation 
I- 

^r^^ - \-^^^ + ^ + t), (1) 

where V/f (r) is the self-consistent Hartrcc potential, 

VH{r^ - mr,t)\^-u / dV' ';: '4 (2) 



The normalization condition for ip is J \ip{r,t)\'^(Pr — N, N is the total number of particles in the condensed state. 
The original Gross-Pitaevskii equation can be obtained by putting u = 0. 

One can write down the Lagrangian density corresponding to this system as follows, 



where * denotes the complex conjugation. The non-linear Schroedinger equation can be obtained from a minimization 
of the action, I = J CSrdt. The BEG of charged Bosons |^ confined in an ion trap can be described by the above 
mentioned Lagrangian if we set —u = e^, where e is the electronic charge. 

To calculate the excitations spectrum of an atomic BEG with gravity-like interaction, we will use the time-dependent 
variational method. This technique has been first used to calculate the low-lying excitations spectrum of a harmonically 
trapped atomic BEG in Ref. The result obtained from the variational method matches with Stringari's result 
within the sum-rule approach. 

In Ref. , it is shown that the oscillation frequencies obtained from the exact ground state and a Gaussian ansatz 
are in good agreement. In this system, a Gaussian ansatz is also a good variational wave function. In order to obtain 
the evolution of the condensate we assume the following variational wave function, 

il^{p,z,t) = C'(t)e-^["W'''+''(*)^'l, (4) 

where C (t) is the normalization constant which can be determined from the normalization condition, p and z are the 
variables in units of A, where A = ^Jh/ mujg = Ti' jmuN is the length scale in this system ( similar to the harmonic 
osscilator length ) when the harmonic trap is absent and lOg ~ mu^N'^ jf? is the gravitational frequency, p is the two 
dimensional vector. a(t) = 1/a^ -t- ia-i and /?(t) = l//3^ -f i(5i are the dimensionless time dependent parameter. ol\ 
and /3i are the condensate widths in a; — j/ plane and along the z direction respectively. The Gaussian variational 
wave function is an exact ground state when the two two-body interatomic interaction is absent. 

Substituting into (^) and integrating the Lagrangian density over the space co-ordinates. We get the following 
Lagrangian, 

^ SNU., ^. 1 2A^ ,1 2 2^ a2a2^ , ^ l^ (1 " )] , , 

i = ^5— [(a?a2 + :j/3?/?2) " (— + o\al) -^(752+ A/^l) - \\-{—r, ^ —)\- (5) 

la I a J ^ Pi V TT a\p\ p\ 



where S = — ""^"''^ is a dimensionless scattering parameter similar to the scattering parameter P — ^ ||l^ for 
an ordinary atomic BEG. Here, oq is the harmonic osscilator length. This S can be positive or negative depending on 
the sign of the scattering length a. The scattering parameter S can also be written as S* = (5287r^/105)/(Afa/AL)^, 
where / = I/Iq and Iq = {4:8n/7){ch'^eQ/maQ)a is the threshold laser intensity to create a self-bound condensate [1). 
For a given intensity / = 1.5, the realistic, stable and self-trapped system ( for sodium atoms ) contains 40 to 400 
atoms 1^ and the corresponding range of 5 is 1 to 100. This range can be alter by changing the scattering length a 

i^[i,l;|;(l — ^)]is the Hypergeomctric fmiction. The last term in the Eq. (||) is the mean- field energy of the 
gravity-like potential. We are interested to find out the excitation spectrum of a self-bound Bose gas as well as in 
the TF-G and G regimes. We have set Vext — because the system is stable even in the absence of an external trap 
potential. 

The energy functional in terms of the variational parameter a in an isotropic system is 
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^ NuS, 3 [2, S 1,, 
By minimizing the energy functional with respect to a, one can get the equihbrium point w which is given by, 



The sound velocity is = fJ^/m, where fi = uS/2a\i/2uP' + 2y^277r(5/zi;^ — l/""^)] is the chemical potential. The 
sound velocity Cs vs the dimensional scattering parameter 5* is shown in Fig.l. 
Using the Euler-Lagrange equation, the time evolution of the widths are. 



«i = 7I + \ -(-3^ + ^^^V^)> (8) 



3 ' \l (9) 



Fa-^ [ai, Pi] is the derivative of F[^, 1; |, (1 — ^)]//3i with respect to ai. Similarly, Fp-^ [ai, /3i] is the derivative of 

1; |, (1 — with respect to /3i. The exact form of Fa^ [ai, /3i] and F/j^ [ai,Pi] are given in the Appendix A. 

We are interested in the low-energy excitations of a gravity-like self-bound Bose condensate. The low-energy 
excitations of the condensate corresponds to the small oscillations of the state around the equilibrium widths ckio and 
/3io- Therefore, we expand around the time dependent variational parameters around the equilibrium widths in the 
following way, ai = aiQ + dai and pi = Piq + 6 Pi. 

The time evolution of the widths around the equilibrium points are 

S S I [2 

.io/3io V TT afo/3fo 2 V TT 



'^"1 ^ + ^Vr 4/3 ~ V 3~752~^^i + o V r^"i["io,/3io,i5ai,(5/3i], (10) 



TT afo/3fo Pto V TT afo/3io V tt 



5Pi = -2\j--^^5ai - (-^ + 2-^/--^^)(5/3i + -Fp,[aiQ, Pw^Sai.Spi]. (11) 



^cei [c^ioj Ao? (^Q!!! <5/3i] is the first order fluctuations around the equilibrium points of Fa-^[ai, Pi]. Similarly, 
F'f3-^[aio, Pio, Sai, SPi] is the first order fluctuations around the equilibrium points of F[}-^[ai, Pi]. The exact form 
of [Q!io, Pio,^o:i, SPi] and i^^^ [aio, Pio, Sai, SPi] are given in the Appendix A. 

We are looking for the solutions of e*"* type. First we solve for these two equation in terms of aio,/3io and later 
we set ckio = Pio = w. For isotropic system, the excitations spectrum are 

3 ^ [2 ^AS 2 



:i-rj + \/z(--— )> (12) 



where Ug = '"'^/^ is the gravitational frequency. The ujj^ and w_ vs. the dimensionless scattering parameter 
S are shown in Fig2. When 5' is small, the gravity-like l/r attractive interaction dominates over the repulsive 
pseudopotential. In this limit, the system becomes more compressible. In other words, the system becomes less 
resistant to density changes. So one would expect that monopole mode lies below quadrupole mode. From the Fig. 2, 
we identify that the upper branch of the excitation spectrum is quadrupole mode ( uj- — ujq ) and lower branch is 
monopole mode ( uj+ = lom )• The monopole and quadrupole modes spectrum in Fig. 2 matches very well with the 
spectrum obtained within the sum rule approach When 5 — —1.169, lum starts decreasing as shown in Fig. 2. At 
Sc = —1.179, the system collapses. The value of Sc can also be obtained from the Eq. (0). For large value of S, 
the monopole mode lies above the quadrupole mode because the repulsive pseudopotential start dominates over the 
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attractive long-range interaction. At S = 17.5, there is a crossing between these two modes which is shown in inset 
of Fig. 2. Interestingly, this crossing of these two modes is also obtained from the time-dependent variational method. 

TF-G regime: For large s-wave scattering length, the kinetic and the trap potential energy can be neglected. 
The gravity-like potential is balanced by the s-wave interaction strength. The total ground state energy is Eq = 

— 0.9648A^^(m/Aj), where Aj = 2Tr\J ah^ /mu is the Jeans wavelength which is the shortest wavelength to keep 
stable condensed state. The ground state energy per particle varies as N. The sound velocity is = fj,/m = 
0.307106(itV5)/m). So the sound velocity Cs varies as N^^^ whereas Cs ^ N^^^ for an ordinary atomic BEC in the 
TF approximation p3| . In this regime, we neglect the contribution of the kinetic energy term in Eqs. (||) and (^) and 
we find that the monopole and quadrupole frequencies are ujm = 0.319951ajg5~'^/^ and ljq = 0.202355(jjgS'^^/^. In this 
regime, the monopole and quadrupole frequencies are obtained by solving the hydrodynamic equations numerically 
in 1^. The monopole and quadrupole frequencies obtained from the variational approach are similar to the exact 
numerical values. For an ordinary atomic BEC in the TF regime, the ujm and ujq are independent of the scattering 
length a |l^Jll[]. But, in this system, it still depends on the scattering length a. Here, the ratio ljm /^q is \/5/2 — 
1.58114. Remarkably, this ratio is identical to the result of Q which is obtained within the sum-rule approach. This 
is also true for trapped atomic BEC without gravity-like interaction in the large N limit. It was first pointed out in 

00. 

G regime: In this regime, gravity-like potential is balanced by the kinetic energy. The trap potential and s-wave 
interaction can be neglected. This is analog of Boson-star (non-relativistic) iQ. The total ground state energy is 
i^o = — (A^/19)?iajg. Using the uncertainty relation, we estimate the total ground state energy is = —{N /\Q)fLLOg 
which is very close to the energy obtained from the variational approach. So our variational ansatz for wave function is 
good in this regime also. The ground state energy per particle varies as N'^. The sound velocity is — 0.159155uS'/rn. 
So the sound velocity Cs varies as iV. In this regime, neglecting the contribution of the s-wave interaction potential in 
Eqs. (H) and (||), we get the monopole and quadrupole modes are ujm = 0.0707355wg and ujq = 0.118363tJg. These 
frequencies are very close to the frequencies obtained within the sum-rule approach j^. Here, ujm/ujq = ■\/5/14 = 
0.597615. This ratio is also identical (up to five decimal) with that of which is obtained within the simr-rule 
approach. 



III. VORTICES OF A GRAVITATIONALLY SELF-BOUND BOSE CONDENSATE STATE 

We consider a gravitationally self-bound Bose condensate state with a vortex along the z-axis. The experimental 
realization of a vortex state would be a direct signature of macroscopic phase coherence of this new atomic BEC with 
a attractive long-range interaction. One can use the time-dependent variational approach to describe the vortex state. 
In the previous section, we have explicitly shown that the monopole and quadrupole mode frequencies obtained by 
using the Gaussian ansatz coincides with the numerical results. So it is natural choice to assume a variational wave 
fimction of a self-bound BEC state with a vortex along the z-axis is, 

ijq{f,t) = Cg{t)p'^e''^^e-'^^^+'^^''>\ (14) 

where q is the vortex quantum number and Cq{t) is the normalization constant. Also, p'^ — + y^, — x'^ + y'^ + 
and 4> — tan^^{y/x). For simplicity, we consider only q — 1 and q — 2. By following the same procedure of the 
previous section, one would obtain the effective Lagrangian which is given by 



(15) 



where gq = J^(q\yi ^ ci = 23/30 and C2 = 37/56. 

The energy functional of the vortex state in terms of the variational parameter a is 

^ NuS ,, 3, , 1 , r2 ,q„S Cn^, , , 

E, = ^[iq+l,){ — ) + J-{^-^)] 16 
Za 2, \ IT a-^ a 

By minimizing the energy with respect to the variational parameter a, one could obtain the equilibrium width Wq 
which is given by 

V2^{q + 1) + j2n{q+lY + l2gqScq 
""^ = 2Z ' 
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where gq and Cq are given above. The system collapses when Sc = —8.53694 for g = 1 and Sc = —25.8875 for q = 2. 
This critical value Sc is increasing with increasing of the number of vorticity. The expectation value of the square 
of the system radius is Iq — v<r^> = N{q + 'i/2)wq. The energy functional satisfies the stability condition, 
s'^^i I ^ n 

The superfluid coherence length, ^, is a distance over which the condensate wave function can heal. In the case 
of a vortex, it corresponds to the distance over which the wave function increases from zero, on the vortex axis, to 
the bulk density. It can be calculated by equating the kinetic energy to the interaction energies. The kinetic energy 
term can not be neglected even for large 5, since it determines the structure of the vortex core. The system exhibits 
superfluid properties if the coherence length is small compared to the size of the condensed state, otherwise it will be 
less prominent to observe the superfluid properties. Equating the kinetic energy to the interaction energies. 



where R — aN^/F / S is the radius of the condensed state at S and, 

i_ s 1 



(18) 



+ V-(— --)• (19) 



is given in the Eq. After rescaling the above equation, we get a quadratic equation of ^ whose solution is, 



e _ f + Vf^ + 8sVf 

R^ SttS ■ 

The superfluid coherence length f/i? vs. S* is shown in Fig. 3 for a wide range of S. The vortex state play an important 
role in characterizing the superfluid properties of Bose system. The critical angular frequency required to produce a 
vortex state is |l5| 



where Eq is the energy of a vortex states with vortex quantum number q and Eq is the energy with no vortex. The 
critical angular frequency fli vs. the dimensionless scattering parameter S is shown in Fig. 4. The critical angular 
frequency decreases with increasing S (or N). For attractive interaction, ilg increases with increasing of S. This is 
also true for an ordinary BEC in the TF regime jlj] . 

The monopole mode frequency for the vortex state with a vortex number q is 

where Wq is given in the Eq. (pT[). 

TF-G regime: For large s-wave scattering length, kinetic energy can be neglected. The superfluid coherence 
length can be obtain from the Eq. (pO|). When S is large, one would gets F = 0.6142S'-i/2 from Eq. (|l|). 
Then, £_/R = 0.17655'-5/«. When is S is very large, coherence length is very small compared to the size of the 
system. The TF-G regime should exhibit the superfluid properties. The mean size of the condensate with q = 1 and 
q = 2 are Ii = 2.21163-\/S' and I2 = 2.4412\/S' respectively. The size of a condensate state with vortices increases 
with the number of vorticity. The critical angular frequencies for q — 1 and q — 2 are Qi = WgO. 00775'^^/^ and 
0,2 = Ci;g0.0094S'~^/^ respectively. The monopole mode frequencies for one and two vortices are uJi — 0.299ujgS^^^^ 
and UJ2 = 0.26u!gS~^^'^ , respectively. These two monopole mode frequencies are less than the monopole mode frequency 
of the vortex free condensate. So, in the TF-G limit, monopole mode frequency of the condensate decreases due to the 
presence of the vortex. The monopole mode frequency for an ordinary atomic BEC in the TF regime is independent 
of the vortex . 

G regime: In this regime, we neglect the contribution from the s-wave interaction energy. The superfluid coherence 
length ^ can be obtained by equating the kinetic energy to the gravity-like interaction energy, /2m^^ = uN/(,, which 
gives ^ = 0.8862_Rg which is almost equal to the radius (Rg) of the condensate state in this regime. As we know that 
if the coherence length is compareable to the radius of the condensate state, it is difficult to exhibit the superfluid 
properties. In this regime, Ii = 12.92 and I2 = 24.84. Here, I2 >> h- The size of a condensate state (with vortices) 
expanding abruptly with increase of number of the vorticity. For example, h/Io ~ 2.8 and I2/I0 ^ 5.39, where 
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Jo = 4.604 is the mean size of the vortex free condensate. The critical angular frequencies for g = 1 and q = 2 are 
ill — 0.0343wg and 0.2 = 0.0215(x'g, respectively. Here, 0,2 < In this regime, the condensate state with vortex of 
q = 2 is unbounded because |/X2|/^ < where /i2 — — 0.0297wS'/a is the chemical potential in the rotating frame. 
So the vortex of q > 2 can not be created in this regime unless there is an additional repulsive potential. Note that 
although there is an indication of instability of vortex with q — 2, this may be just an artifact of the variational 
approach. The monopole mode frequency for q = \ is uji = Q.QlAQujg. The uoi is also less than the monopole mode 
frequency lom in the vortex free condensate. 

It should be noted that the vortex has two different length scales, condensate radius and core radius, whereas the 
trial wave function (14) has only one variational length scale (a). In this variational approach, the various numerical 
values computed for the energies of the two regimes and for the collapse are just indicative but more accurate values 
can be obtained by other rigorous methods. 



IV. SUMMARY AND CONCLUSIONS 



In this paper, we have derived an analytic expressions of the monopole and quadrupole excitation frequencies of a 
self-bound Bose gas induced by the electromagnetic wave for a wide range of the dimensionless scattering parameter 
S. Later, we consider the two new regimes, namely, TF-G and G regimes. In these regimes, we have calculated the 
lower bound of the ground state energy, sound velocity, monopole and quadrupole mode frequencies. Our results are 
in excellent agreement with the result obtained by using the sum-rule approach. Interestingly, the ratio lom I^q is 
identical (in both the regimes) with that of which is obtained within the sum-rule approach. In the TF regime of 
an ordinary atomic BEC, the monopole and quadrupole mode frequencies are independent of the scattering length a. 
On the other hand, in the TF-G regime, the monopole and quadrupole mode frequencies depends on the scattering 
length a. The local sound velocity Cs varies as 

in the TF-G regime, whereas Cs ~ N^^^ for an ordinary atomic 
BEG in the TF regime. For harmonic trapped Bose system, the excitation frequencies are determined by the oscillator 
frequency of the trap potential. But, in this system, the monopole and quadrupole mode frequencies are fixed by the 
gravitational frequency lOg. 

In section III, based on the time-dependent variational method and simple ansatz for the wave function (equ. [ij ), 
we have studied a rotating gravity-like self bound BEG states with vortices along the z-axis. We derived an analytic 
expressions for the coherence length and the critical angular frequency to create a vortex in the condensed state. We 
found that the coherence length in the TF-G regime is very very small compared to the radius of the system. On the 
other hand, the coherence length in the G-regime is compareable to the radius of the system. We could say that the 
TF-G regime should exhibit the superfluid properties more prominently than the G-regime. In the TF regime of an 
ordinary atomic BEG, the monopole mode frequency of the condensate does not change due to the presence of the 
vortex. But, the monopole mode frequency in the TF-G regime as well as in the G regime of this new BEG decreases 
due to the presence of the vortex. 

I would like to thank G. Baskaran and Pijush K. Ghosh for valuable comments. 



APPENDIX A: 



Here, we shall give the exact form of F[ai,/3i] = 1; |; (1 — ^)]//3i, its derivative with respect to ai, /3i and 
the first order deviation from the equilibrium functions [aioj /3io] f^nd Fp^ [aio: /3io]- 



The derivative of i^[Q;i,/3i] with respect a\ is given by 

Expanding Fa^[ai, Pi] around the equilibrium width aio and Piq. The first order deviation from the equilibrium 
bmction -Fqi [aio, /3io] is. 



F^, [aio, Pio, Sai,SPi] = ^[-|^ + i(-^ + 12^) + ^(-J- + 36^ - 30^) + ]6ai (A3) 

Pio 3Pio 5 Pio Pio ' Pio Pio Pio 
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1 , 2 aTn 1, Q!in a?n N 1. Q^in Ctfn "f n ^ 

HlO A-'lO ^ A-'lO f-^lO ' A-'lO yiO PlO 

PIO o Pio Pio ' Pio Pio Pio 



The derivative of F[q;i,/3i] with respect to /3i is given by 



„ r „ 1 r 1 / a? , 1 , a? af , 1 , a? af a? , , , . 

^^^f"- = -^[^ + 3^^ - ]f ) + 5^^ ~ + 7^^ - 3^ + 3^ - ^) + ] (A4) 

Expanding Fp-^^ [ai, Pi] around the aio and /3io. The first order deviation from the equihbrium function Fp-^^ [aio, Pio] 
is the foUowing, 

i^0jaio,/?io,foi,<5/3i] = --^[---H + (_4^ +4 10) + ]Sai (A5) 

Pio Pio ^ Pio Pio 



PlO O Pio Pio Pio ' Pio Pio Pio 



73r ii + o(i - + ■^^^ " + 7^) + 7(1 - ^^7^ + -^7^ - 7^) + l^Pi 

PlO 'J PlO ^ PlO PlO ' PlO PlO PlO 

I + 1-5/31 

20|^^) + i(-18^ + 60|^l-42|^I) + 

PlO ' PlO PlO PlO 
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FIG. 1. The sound velocity Cs as a function of the dimensionless scattering parameter S. 
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FIG. 2. Monopole and quadrupole mode frequencies vs. the scattering parameter S, for positive as well as negative values 
of a. The solid and dashed lines corresponds to the quadrupole and monopole mode frequencies, respectively. The crossing 
between two modes is shown in the inset of this figure. 
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FIG. 3. The superfluid coherence length ^ as a function of the dimensionless scattering parameter S. 
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FIG. 4. The critical angular frequency Qi as a function of the dimensionless scattering parameter 5*. 



